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By A. Trouve 
!• Motivation and objectives 

In the flamelet theory, the modeling of the mean reaction rate is based primarily 
on a statistical description of the wrinkling of the flame surface due to the turbulent 
motions. The amount of wrinkling is quantified by the mean flame surface- to- 
volume ratio, also called the flame surface density. The flamelet theory produces 
an exact evolution equation for the flame surface density, called the E-equation, 
where basic physical mechanisms such as production by hydrodynamic straining and 
dissipation by propagation effects are described explicitly. In flamelet models, the 
flame surface density is obtained via a modeled formulation of this exact E-equation 
(Marble & Broadwell 1977; Darabiha e* at. 1987; Pope & Cheng 1988; Maistret et 
at. 1989; Candel et at. 1990; Cant, Pope & Bray 1990; Borghi 1990; Mantel & Borghi 
1991). Recent work at CTR, aimed at improving current models for the E-equation, 
focused primarily on the production and dissipation terms due to the turbulent flow 
field (Meneveau Poinsot 1991; Trouve & Poinsot 1993). In these studies, the 
turbulence upstream of the flame is considered isotropic and the mean flow velocity 
gradient tensor, \dUi / dx j ] , features only one non-zero component that corresponds 
to the mean flow acceleration through the combustion zone. Flame configurations 
featuring more complex mean deformation patterns, as found in practical systems 
where the combustion occurs in shear layers, boundary layers, or stagnation point 
flows, have not been considered yet. 

There is, in fact, little knowledge on the influence of mean velocity gradients on 
the dynamics of premixed flame surfaces. In current flamelet models, since the tur- 
bulent flame stretch in the equation for E depends on the turbulent kinetic energy, 
&, and its rate of dissipation, 6, the effects of mean deformations are incorporated 
implicitly through their impact on k and e. An additional term that is linearly re- 
lated to [dUi/dxj] is also sometimes introduced to account for possible direct effects 
(see Duclos, Poinsot & Veynante 1993). 

In the present work, we use three-dimensional Direct Numerical Simulation (DNS) 
of premixed flames in turbulent shear flow to characterize the effect of a mean shear 
motion on flame surface production. The shear is uniform in the unburnt gas, and 
simulations are performed for different values of the mean shear rate, 5. The data 
base is then used to estimate and compare the different terms appearing in the 
S-equation as a function of S. The analysis gives in particular the relative weights 
of the turbulent flow and mean flow components, a? and A T , of the flame surface 
production term. This comparison indicates whether the dominant effects of a mean 
flow velocity gradient on flame surface area are implicit and scale with the modified 
turbulent flow parameters, k and e, or explicit and scale directly with the rate of 
deformation. 
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1.1 The evolution equation for the flame surface density 
In the flamelet regime, it is convenient to describe the flame-flow interactions in 
terms of the following two basic ingredients: a flame speed that characterizes the 
mean combustion intensity, and a surface area that characterizes the flame won mg 
(Williams 1985; Peters 1986; Bray 1990; Pope 1990; Poinsot, Veynante & Can del 
1991). For instance, the mean reaction rate may be written as a mean integrate 
chemical rate times the flame surface density: 


(wfl) = ( J a>Rdn)s (S') — ( J WRdn)s S, 


( 1 ) 


where ur is the mass of fuel consumed per unit time and per unit volume; / URdn 
is the local integral of the reaction rate along the flame normal direction; E is 
the flame surface area per unit volume; E is the flame surface density, define as 
the expected value for E\ E = (S'). Note that in (1), all quantities are local and 
mean operators correspond to ensemble-averaging. While ( ) denotes a standard, 
unweighted ensemble-average, ( )s denotes a surface mean, define “ an area " 

weighted ensemble-average (Pope 1988), (Q)s = (Q s )/(£)- (Q S )/ S - 

In (1), the integral of the reaction rate can be replaced by the local fuel consump- 
tion speed, 5c, defined as: 


5c 


Pu^ R,u J 


u>ndn , 


( 2 ) 


where p u and Yr iU are respectively the density and the fuel mass fraction in the 
unburnt gas. The mean reaction rate may then be expressed as the surface mean 
of 5c, called the flamelet speed, times the flame surface density: 


(w«) = (/> u Vr, u (5c)s) £. 


(3) 


(3) is the classical flamelet expression for the ensemble mean of the reaction 
rate. In this expression, the flamelet speed, (5c)s, accounts for local variations of 
the reaction rate along the flame surface. Recent direct simulations suggest that 
for flame Lewis numbers close to unity, the departures of ( S c )s from the laminar 
consumption rate tend to average out when integrated across the whole turbulent 
flame (Ashurst, Peters k Smooke 1987; Haworth k Poinsot 1992; Rutland k Trouve 
1993; Trouve k Poinsot 1993). In this situation, the flame surface density is to first 
order the single key quantity that determines the mean reaction rate. Note that 
Trouve k Poinsot (1993) also describe one flame with thermal-diffusive instability 
where the effect of the turbulence on the flame surface area is coupled with a 
significant increase in the mean flame intensity. In that case, both the flamelet 
speed and the flame surface density need to be modeled to determine the mean 

reaction rate. . n . 

In many situations, the principal effect of the turbulence is for the fluctuating 

velocity field to wrinkle the flame and greatly increase its surface area. This phe- 
nomenon is represented in (3) by the flame surface density, E. Following the pioneer 
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work of Marble & Broadwell (1997), an exact balance equation for the flame surface 
density was derived by Pope (1988) and Candel & Poinsot (1990): 

as 

+ ^-(X) 5 S = (k) s E, (4) 

where X is the velocity of the flame surface, given by the sum of the fluid velocity 
and the flame propagation speed in the normal direction: X = u + um; n is the unit 
vector normal to the flame surface; and k is the flame stretch. A useful expression 
for k is in terms of strain rate, flame curvature, and flame propagation speed (Candel 
& Poinsot 1990): 

k = a T + 2wk m , (5) 

where ay is the rate of strain acting in the flame tangent plane: ar = V.u-nn : Vu 
(we use here standard tensorial notations: nn : Vu = n.njdu./dxj)-, and k m is the 
flame surface curvature, as given by the divergence of the flame normal direction: 
2 km = V.n. In (5), positive curvature is chosen convex towards the reactants. 

Equation (4) can be cast in various forms. For modeling purposes, it is useful 
to split the velocity vector into a mean component and a turbulent fluctuation: 
u = U + u", where the tilde superscript denotes a Favre-averaged quantity: U = 
(pU)/(p)- We can then re-write (4) as follows: 

as 

~fa + V.US + V.(u %E + V.(u>n) s E = (. A T )s E + (a T ) s S + 2 (wk m ) s E, (6) 
where we use the following notations: 

(A t )s = VU - (nn) s : VU, (7) 

(a T )s = (V.u" - nn : Vu") s . (8) 

The three convective terms on the left-hand side of (6) are transport terms that 
correspond respectively to convection by the mean flow, turbulent diffusion, and 
flame propagation. The terms on the right-hand side of the equation are the source 
and sink terms for the flame surface density: ( At)s is the tangential strain rate 
due to the mean flow field; ( a T )s is the strain rate due to the turbulent flow field; 
and 2 {wk m }s is a term that accounts for the combined effects of flame curvature 
and flame propagation. 

IVe are particularly interested in this study in the effect of mean flow deformations 
on the dynamics of flame surface densities.^ As seen in the equation for E, in the 
presence of mean flow velocity gradients, dUi/dxj, the flame surface is modified in 
two different ways: first by a rapid distortion that accounts for the direct straining 
of the flame surface by the mean velocity gradients (this effect is represented in (6) 
by ( A t )s ) and second by a slower modification resulting from the adjustment of the 
turbulence to the applied deformations and from the corresponding variations in the 
turbulent flame stretch (this effect is represented mainly by ( ax)s and 2(wk m )s). 
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1.2 Direct numerical simulation of turbulent premixed flames 
We use DNS to analyze the different terms in the equation for the flame surface 
density. The simulations are performed using a three-dimensional, compressible 
Navier-Stokes solver that fully resolves the turbulent flow field. The solver uses a 
high-order finite difference scheme: spatial derivatives are computed with a modified 
Pade scheme that is sixth-order accurate (Lele 1992); solutions are advanced in 
time with a third-order Runge-Kutta method (Wray 1990); boundary conditions 
are specified with the NSCBC method (Poinsot k Lele 1992). We refer the reader 
to Lele (1992) and Poinsot k Lele (1992) for further details concerning the system 
of equations solved and the numerical methods. . , 

Because of the otherwise prohibitive computational cost, simulations are limited 
to simple but finite-rate reaction schemes. In this work, the chemistry model is a 
single step, irreversible chemical reaction where the reaction rate depends exponen- 
tially on temperature (Arrhenius kinetics): 

wr = BpYr exp ^ ^ 

where T„ is the activation temperature and B is a constant that is determined 
according to the selected laminar flame speed. This formulation corresponds to a 
binary reaction in which one of the reactants, Yr, is strongly deficient, for example, 
in fuel-lean combustion. Also, it is worth emphasizing that the simulations do not 
use the constant density assumption; the reaction is exothermic and heat release 

effects are fully accounted for. 

Following Williams (1985), we re-write the reaction rate as: 

/ -/3(l-en nm 

= ApYR exp 

where 0 is the reduced temperature, 0 = (T-T u )/(T t -T„); T u is the temperature 
of the fresh reactants; T b is the adiabatic flame temperature; and the coefficients 
A, a, and /? are, respectively, the reduced pre-exponential factor, the heat release 
factor, and the reduced activation energy: 

A = Bexp(— /3/a), c* = (T b — T u )/Tb, and fi = aT a /Tt,. (11) 

The values of the flame parameters are reported in Table 1. 

Another important feature of the simulations is that the transport coefficients 
are temperature dependent. These coefficients satisfy the following relations: 

p = p u (T/T u ) h , Le = \th/pDc p = constant , Pr = pc p /\,h = constant , (12) 

where p, \ t h, and D are the molecular diffusivities of, respectively, momentum, 
internal energy, and species mass; c p is the specific heat at constant pressure; b 
is a constant; and Le and Pr are respectively the Lewis number and the Prandtl 
number. The simulations were performed for a non-unity Lewis number, Le = 0.8. 
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Figure 



Le a 0 A+ b Pr 


0.8 0.75 8.0 456 0.76 0.75 


TABLE 1. Dimensionless flame parameters, A is made non-dimensional 
by the laminar flame time, A + = A (A th/pc p ) u /sL 2 . 

The computational configuration corresponds to a premixed flame embedded in 
a three-dimensional, turbulent shear flow (see figure 1). The calculations are ini- 
tialized with fresh reactants on the top side of the domain (x 2 > 0) and burnt 
products on the bottom (x 2 < 0); the two are separated by a plane laminar flame. 
The turbulence is initially isotropic, its velocity field being specified according to a 
model energy spectrum: 

m = c e* P f-2(i) ! >, (is) 

where k is the wavenumber and C and k t are model parameters that are specified 
according to, respectively, the initial turbulent rms velocity, u", and integral length 
scale, /. 

At t = 0, a simple shearing motion, 0U\ j dx 2 = S, is applied on the flow system. 
Note that due to the flow acceleration through the flame, the shear rate does not 
remain uniform in the cross-stream x 2 direction; in the flame zone, the velocity 
gradient dU i / dx 2 is a function of both time and cross-stream position. Outside the 
flame zone, however, the turbulence remains homogeneous in all three directions. 
In the burnt gas, the mean velocity gradient changes with time; it remains constant 
and equal to S in the unburnt ga s. 

The top and bottom walls of the computational domain are inflow and outflow 
boundaries while periodic boundary conditions are applied in the x, and x 3 direc- 
tions. Non-periodic boundary conditions are needed along x 2 for the following two 
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reasons: to be able to maintain numerically a constant shear; to account for the 
density change from unburnt to burnt gases. Note, however, that no turbulence 
is generated at the inflow boundary, and the simulations are time-evolvmg rather 

than space-evolving. . , , 

Values of the run parameters are reported in Table 2. At time t - 0, the tur- 
bulence is characterized by flow length scales that are a few times larger than the 
laminar flame thickness, l/l F » 4 - 5 with l F = (W pc p )»/s L , and turbulence 
intensities that are higher than the laminar flame speed, u"/s L > 5 . The initial 
influence of the mean shear on the overall turbulence dynamics is measured by 
a time scale ratio, 5r 0 , where r 0 designates the initial, turbulent eddy turn-over 
time, ro = l/u". All simulated cases correspond to similar values of the miti 
Reynolds number. The initial turbulent Reynolds number based on the Taylor mi- 
croscale, A, is approximately Re\ « 25. The initial turbulent Reynolds number 
based on the integral length scale is Re, « 35. In Table 2, a Damkohler number 
is also introduced using the Kolmogorov time scale and the laminar flame time: 
Da = (A /u")/(l F /s L ) = ViE(Ti k /v k )/(l F /sL), where r] k and v k are respectively the 
Kolmogorov length and velocity scales. While the simulated flames correspond to 
relatively low Damkohler numbers, Da < 1, the reaction zone can still be described 
as a surface and the simulated flames occur in the flamelet regime. 

Since in the presence of mean shear the turbulent structures are elongated in the 
streamwise direction, we use a computational domain that is twice as long in t e 
streamwise direction compared to the cross-stream and spanwise directions. The 
grid spacing is uniform in all three directions; the grid resolution is 129 x 129 x 65 
and is, therefore, twice better in the direction of mean flame propagation, x 2 ■ Also, 
given the growth of the turbulent length scales in the streamwise direction, the 
numerical simulations are all limited to relatively low total shears, St < 15, after 
which the eddies have sizes that are comparable to the computational domain and 
the numerical simulations suffer from insufficient resolution. 


Run designation 

St o 

Re, 

Re\ 

Sl F /si 

vT/sl 

l/l F 

Da 

Le 

Case 1 (lew08tSR) 

10.5 

35. 

25. 

10. 

5.0 

5.3 

0.8 

0.8 

Case 2 (lew08SR) 

7.0 

35. 

25. 

6.7 

5.0 

5.3 

0.8 

0.8 

Case 3 (lew08sSR) 

3.5 

35. 

25. 

3.3 

5.0 

5.3 

0.8 

0.8 
A O 

Case 4 (lew08wSR) 

3.3 

38. 

28. 

6.7 

7.5 

3.8 

0.4 

0.8 
A A 

Case 5 (lew08uSR) 

1.7 

38. 

28. 

3.3 

7.5 

3.8 

0.4 

0.8 

A A 

Case 6 (lew08vSR) 

0.7 

38. 

28. 

1.3 

7.5 

3.8 

0.4 

0.8 


TABLE 2. Dimensionless parameters for the reacting flow simulations. 

All parameters correspond to the initial condition, t = 0. 

All terms in (6) can be obtained from the simulations. We refer the reader 
to Trouve (1993) and Trouve & Poinsot (1993) for further details concerning the 
diagnostics that were developed to extract the relevant information from the DNS 
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Figure 2. 

the relative 
2 and 6. 



database. We simply briefly recall here that flame-based quantities are defined 
by making the assumption that the resolved finite-thickness reaction zone can be 
viewed as a constant progress variable surface, c = c f = 0.8wherec = 1-(YrIY r<u ). 
Ensemble- averaging is conveniently performed in the homogeneous ii - r 3 planes. 
For completion, we recall the expressions used to calculate the flame propagation 
speed, the flame surface density, and the surface means: 


w = X.n — u.n = 


1 Dc 


E = 


|Vc| Dt p\vy r \^ ur v (/ , - dv *r))> 

- — / 

Jc=Cf 


|Vc| 




dl, 


(14) 

(15) 


where L Xl and L X3 are the x i and x 3 dimensions of the computational domain. 


(Q)s = 


(gg) 

(E') 


+ (&)*) 

(l v 'l ■«/\/<£r) 2 + (&>*) ' 


(16) 


2. Accomplishments 

The simulations describe the wrinkling of the flame zone by the flow as well 
as the transition of the flow field from isotropic decaying to anisotopic sheared 
turbulence. The mean shear rate, S , is used as a control parameter to modify the 
coupling between the mean flow, the turbulence, and the flame. Since the wrinkling 
of the flame is mainly driven by the turbulent motions in the unburnt gas, we first 
describe in §2.1 the mean shear- turbulence interactions that modify the turbulence 
upstream of the flame. We then turn in §2.2 to the impact of these modifications 
on the evolution of flame surface densities with emphasis on the effect of S on the 
flame surface production term. 
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2.1 The effect of mean shear in homogeneous turbulent shear flow 

Turbulent shear flows present several additional levels of complication compared 
to isotropic turbulence because of the anisotropy of the flow, the production of 
turbulent kinetic energy by the mean strain, the possible amplification of hydro- 
dynamic (Kelvin- Helmholtz) instabilities, and the resulting presence of large scale 
coherent structures. The general complexity of sheared turbulence is significantly 
reduced in the case of homogeneous shear flows, where the mean shear is uniform 
throughout the flow. For instance, such basic flow configurations are free of large 
scale hydrodynamic instabilities. In addition, homogeneous turbulent flows allow 
much simplified and more powerful statistical descriptions. Due to this increase m 
simplicity, homogeneous shear flows have been extensively studied in the literature, 
experimentally (Tavoularis & Corrsin 1981; Rohr et al. 1988; Tavoulans & Karnik 
1989), as well as numerically (Rogallo 1981; Rogers & Moin 1987; Lee, Kim & Mom 
1990; Holt, Koseff & Ferziger 1992). 

We are particularly interested in this work in the effect of changing the mean shear 
rate, 5, on the energy levels and the structure of the turbulence. To characterize the 
effect of 5, we performed a preliminary series of six simulations without flame, lhe 
run parameters are identical to the ones given in Table 2, except for the absence of 
chemical reaction. The grid resolution is 129 x 65 x 65. The results are used m the 
subsequent simulations with flame to describe the dynamics of the turbulence in the 
unburnt gas. We present in this section the main results from these non-reacting 


flow simulations. 

The evolution of initially isotropic turbulence in a mean flow of uniform sim- 
ple shearing motion, 5 = dUJdx 2 = constant , can be described in terms of the 
following three dimensionless parameters: the total sheax, St , that gives a non- 
dimensional time; the time scale ratio 5r 0 ; and a turbulent Reynolds number, 
Rei = u"l/v, taken at t = 0. 5r 0 is a measure of the initial importance of the 
applied shear. Re t determines the impact of viscous effects. Both parameters serve 
to describe the effect of initial conditions, and their influence will be negligible at 


large times. 

The present simulations, however, are limited to total shears of approximately 
10. At time St - 10, the turbulent eddies remains in non-equilibrium. For instance, 
the quantity Sk/e is still evolving. It is generally believed that uniformly sheared 
turbulence will ultimately achieve a self-preserving state in which the turbulence 
changes at a rate independent of the initial conditions. In the simulations, the 
turbulence is clearly still far from this equilibrium state. 

For small values of the total shear St, the simulations reveal a range of flow 
regimes ranging from viscous decay for small values of Sto, to rapid distortion for 
St 0 > 5. This diversity in the turbulent flow response to the mean shear is shown in 
figure 2. While for St 0 = 0.7 the flow goes through an initial phase characteristic of 
decaying grid turbulence and energy growth is only observed at later times, 5r 0 > 3, 
for St 0 = 7, the flow is dominated by the mean shear-turbulence interactions and 
the energy growth starts immediately at t = 0. 

Figures 3 and 4 describe the partition of turbulent kinetic energy between the 



237 


Premixed flames in turbulent shear flow 



total shear ; St 


Figure 3(a ). Pa rtition of the turbulent kinetic energy, k\ if,, = u n u n /2k* 
Q 12 = ^ 12 / \ZSTii A^ 22 - Case 6, 5ro = 0.7. * J 



Figure 3(B) . Partition of the turbulent kinetic energy Jfc- Jb, 
Case 6, 5r 0 = 0.7. ’ 


«/2 = x *. 


T d S P“ wis ' Veiocit >’ As seen in flgures 3(a) 

( ), simulations feature an initial transient phase, 0 < St < 6 where the 

SiuslZnttr 18 ’ Chan « e ra P idl y> Allowed * a phase of slower 

justment where the flow approaches structural equilibrium with Ku asymptoting 

Wnstaj, values Note that the flow remains more isotropic in case 6 coLp^d "o 

sWt'JAn T con ! ls ‘ ,!nt with *>>e findings by Lee, Kim k Moin (1990), who 

wh™ the ,^ 7 f Sh < f"' driven fnrbulence changes dr asically 

flo “ir^s w (,h “ ’ * “ mCre “ ed “ d ‘ hat “» d ^“ * «*— <* tta 

Figures 3(b) and 4(b) compare the magnitudes of the different turbulent kinetic 
energy components k, = u ''u'> / 2 = K it x k. The initial decay phase is very pro- 
nmmced m c^se 6t m flgur. 3(b), the different turbulent inteLities remjn lower 
than then initial value up to S. = 8, The situation is drastically different at high 
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Figure 4(a). 
Q 12 = Kn/y/Rn 


Partition of the turbulent kinetic energy, k; K,j 
A' 22 - Case 2, Stq = 7. 


*>y/2Jb; 



FIGURE 4(b). Partition of the turbulent kinetic energy, t, fc, 
Case 2, 5r 0 = 7. 


u"u" j2 = K« x fc. 


shear rates: in figure 4(b), at St = 8, the streamw 1S e component h is already an or 
der of magnitude larger than its initial value. Furthermore, the simulations feature 
the following classical arrangement: ki > fa > k 2 . The first inequali y re ec s 
fact that the streamwise component receives energy directly from the mean shear, 
while the other components are maintained by redistribution of energy through the 

pressure-strain correlation. , i . u 

We also present in figure 4(b) a comparison between the DNS data and the pre- 
dictions from Rapid Distortion Theory (RDT) (see for example Townsend 1976; 
Maxey 1982). RDT is a linear theory that applies to situations in which the mean 
shear-turbulence interactions dominate the nonlinear turbulence- turbulence inter- 
actions. Its domain of validity corresponds to 5r 0 > 1. RDT sheds light on e 
mechanisms of energy exchange between the different turbulent velocity compo- 
nents. For instance, RDT shows that while most of the energy is concentrated in 
the streamwise component, k u energy is also transferred from both fci and k 2 to k 3 
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by the pressure-strain correlation, which accounts for the net reduction (increase) 
in &2 (£3) by the shear distortion seen in figure 4(b). 

Hence, while the initial time evolutions of ki and k$ depend strongly on the 
parameter Stq, the cross-stream turbulent velocities are less sensitive to its value: 
for Stq < 1, A?2 decreases at small times because of viscous decay; for Stq 1, jt 2 
also decreases because of a net energy loss to /fc 3 . 

2,2 The evolution of flame surface area in sheared turbulence 

It is well known that in the presence of mean shear, the turbulent rates of mixing 
and chemical reaction tend to be increased. One reason for this increase is the devel- 
opment of large scale coherent structures that are amplified by the Kelvin- Helmholtz 
hydrodynamic instability. This pseudo-laminar mechanism has been extensively 
studied in the last twenty years, including in premixed and nonpremixed combus- 
tion systems. In homogeneous turbulent shear layers, the hydrodynamic instability 
does not occur, and the enhancement of mixing and chemical reaction results from 
increased straining and increased turbulent transport due to shear-production of 
turbulent kinetic energy. This turbulent mechanism has been studied in a number 
of recent numerical studies in the context of mixing problems (Rogers, Mansour 
Reynolds 1989; Holt, Koseff & Ferziger 1992) and nonpremixed chemical reactions 
(Leonard Hill 1992; Nomura & Elghobashi 1992). Nomura & Elghobashi assume 
in their study an infinitely fast reaction, while Leonard & Hill consider a finite-rate 
reaction both with temperature-independent and temperature- dependent kinetics. 
Both studies are limited by the constant density assumption. In these studies, the 
enhancement of chemical reaction due to the mean shear is explained by the in- 
creased turbulent strain rates and increased reactant concentration gradients found 
along the flame surface. Since in the absence of quenching, higher concentration 
gradients are equivalent to higher reaction rates, it is suggested that one basic effect 
of mean shear is to increase the mean flame intensity. This effect is related in (3) to 
modifications of the flamelet speed, (Sc) s, not to modifications of the flame surface 
density, E. 

In premixed flames, however, the departures of (Sc)s from the laminar consump- 
tion rate tend to average out when integrated across the whole turbulent flame, and 
in many situations the mean combustion intensity is only weakly sensitive to the tur- 
bulence. Therefore, the argument that relates the increased combustion efficiency 
observed in nonpremixed flames with shear to a change in the flame structure is not 
expected to apply to premixed configurations. In that case, the effect of mean shear 
must be related to a modification of the total flame surface area. In the following, 
we use DNS to get basic information on the dynamics of flame surface densities 
as a function of the mean shear rate, S. The mean shear rate, S', is varied in the 
proportion of 1 to 10, as seen in Table 2. Large values of the shear rate (as in 
cases 1, 2, and 4) are relevant to a flame propagating in the near- wall region of a 
turbulent channel flow. Lower values (as in case 6) can be found in flames stabilized 
in subsonic shear layers. 
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Figure 5(a). The effect of the turbulence-to-mean shear time scale ratio, 5r 0 , 
on the relative increase of total flame surface area, ( Sv)/L Xl L Xa . A comparison 
between cases 4, 5, and 6. Time is made non-dimensional by the initial, turbulent 
eddy turn-over time, tq. 



total shear, St 

FIGURE 5(b). The effect of the turbulence-to-mean shear time scale ratio, St 0 , 
on the relative increase of total flame surface area, ( Sv)/L Xl L X3 . A comparison 
between cases 4, 5, and 6. Time is made non-dimensional by the mean shear rate, 

5. 

2.2.1 The total flame surface area 

As described in §1.2, six different cases have been simulated that correspond to 
turbulent flames characterized by the same laminar flame thickness, If , the same 
laminar flame speed, s L , embedded in two different initial turbulent flow fields (in 
this respect cases 4, 5, and 6 differ from cases 1, 2, and 3), and with different values 
of the mean shear rate, S. 

The effect of changing the shear rate on the production of flame surface area, 
(S v ), is shown in figure 5 for cases 4, 5, and 6. In figure 5(a), (S v ) is plotted 
versus time, with time made non dimensional by To. Cases 4 and 5 feature an 
initial phase of slow growth followed by a phase of more rapid growth. In case 6, 
the transition to this second phase is not observed. These two successive phases are 
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more clearly seen when (S v ) is plotted as a function of total shear, St (figure 5(b)). 
For 0 < St < 10, the production of flame surface area is slow and the differences 
between all simulated cases remain moderate. In contrast, for St > 10, the flame 
surface is rapidly growing at a rate that does not seem to depend on S. In case 6, 
the simulation is stopped at St ~ 8.5, before the transition to the second phase is 
completed. 

One important result in figure 5(b) is that the flame response to the applied shear 
scales like the turbulent flow field (§2.1) with the total shear, St. 

2.2.2 The production and dissipation of flame surface area 

Using the diagnostics presented in §1.2, we now turn to the analysis of the source 
and sink terms of flame surface area as they appear on the right-hand-side of equa- 
tion (6). These terms can be resolved spatially across the turbulent flame (TVouve 
1993, Trouve & Poinsot 1993). In the following, they are space-averaged in the 
cross-stream direction, and we limit our discussion to time variations of global 
flame properties. We estimate the global value of any quantity Q using the follow- 
ing E-weighted space-averaging scheme: 

= ( J (Q)s(x 2 ,t) E (i 2 ,f)dr 2 )/( J ^{x 2 ,t)dx 2 ). (17) 

Using^(17), we define the total production of flame surface by hydrodynamic stain* 

P i total dissipation by flame propagation effects, D y and the resulting net 
total flame stretch, k as follows: 


P(t) = At + clt, 


(18) 


D(t) = 2wk m , 
K(t) = P + D. 


(19) 

( 20 ) 


Also, the total flame surface area, (5V) 5 is equal to the volume- integral of the 
flame surface density: 

^ = J Z(x 2 ,t)dx 2 , (21) 


and, using (4), k can be directly related to the instantaneous rate of change of the 
total flame surface area in the computational domain: 


= ( 22 ) 

Figure 6 shows P, D, and k as a function of total shear for cases 4, 5, and 6. It 
is seen that while the simulated cases exhibit large differences in the time history 
of the production term, P, these differences tend to be balanced by corresponding 
variations in the dissipation term, D, and the net total flame stretch, k, follows 
approximately the same time evolution for all cases. In other words, the total 
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FIGURE 6. Time evolution of flame surface production, P, flame surface dissipa- 
tion, D, and total flame stretch, k. A comparison between cases 4, 5, and 6. P, D, 
k are made non-dimensional by the laminar flame time, (A** / pCp)u/ s L ■ 



FIGURE 7. ot/At vs St. Relative weights of ar and At in their contributions to 
flame surface production, P. A comparison between cases 4, 5, and 6. 

flame stretch scales with the total shear, St, but does not depend on the parameter 
5r 0 . It is worth emphasizing that this result is somewhat unexpected given the 
large impact of St 0 on the turbulence dynamics as described in §2.1. One quantity, 
however, that remains approximately the same from one simulation to the other is 
the cross-stream turbulent rms velocity, as measured by k 2 . It is believed that the 
insensitivity of the total flame stretch to changes in Stq is a strong indication that 
k scales with rather than k. 

Furthermore, since k appears to be a function of St and not of Sr 0n (22) can be 
re-written as: 

(Sv)(St) = ( S v )(t = 0) exp k (T) dT ^ . (23) 

(23) shows that at constant total shear, St = constant, the flame surface area is a 
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0 5 10 15 20 

total shear, St 


FIGURE 8(a). The effect of S on the magnitude of — nin2. A comparison between 
cases 1 to 6. 



total shear, St 

FIGURE 8(b). The effect of S on the magnitude of At 12 - A comparison between 
cases 1 to 6. At 12 is made non-dimensional by the mean shear rate 5. 

decreasing function of 5, as was found in figure 5(b). 

As seen in (18), the mean tangential strain rate at the flame can be decomposed 
into two components: At represents the contribution of the mean rate of deforma- 
tion; and clt represents the effects of the turbulent flow velocity gradients. Figure 
7 presents the ratio ot/At versus non-dimensional time for cases 4, 5, and 6. At 
0, the flame is plane and there is no strain rate due to the mean flow field, 
At = 0. As thejflame gets wrinkled, the effect of the mean flow field gets more 
pronounced and At and ar take comparable values. For St > 5, however, both the 
absolute value of ar and its relative weight in the expression for P are increasing. 
These results indicate that a transition to a fully turbulent regime is occurring: 
for St > 5, the turbulence levels, k\ , £3, and also k2 are all increasing; the hy- 
drodynamic straining at the flame, as measured by P, takes higher values; in the 
expression for P, clt is increasing faster than At and the turbulent flow component 
of flame production is becoming dominant. 
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We now take a closer look at the strain rate due to the mean flow field. (At)s 
is a production term that depends on the mean cross-stream velocity gradient and 
the flamelets orientation: 

(A T )S = (1 - (”2 2 )s) ~ ( nin2 )s^> ( 24 ) 

where n, is the x, -component of the flame normal vector. The first term in (24) 
accounts for the flow acceleration through the flame; the second term describes 
the rapid distortion of the flame surface due to the applied shear. Both terms 
remain positive across the turbulent flame. Note that since in (24) the velocity 
gradients are multiplied by quantities that are direct measures of the turbulent 
flame wrinkling, the magnitude and relative weights of the two components of ( At)s 
are difficult to assess a priori . In the simulations, the statistics describing the 
flamelet’s orientation do not vary significantly across the turbulent brush; (n 2 2 ) s 
and (nin 2 ) s are functions of time but nearly independent of x 2 position. Therefore, 
space-averaging across the turbulent flame yields the following approximation: 


At = At 22 + At 12 


dU 2 ^ „dlh 

“&7 + /, aV 


(25) 


where a and j3 are constant coefficients. 

Figure 8(a) presents the time variations of -nTn 2 which can be considered as an 
estimate of /3. All simulated cases follow the same time evolution when scaled in 
terms of St: growth from the initial zero value; saturation at St ~ 5; short period 
of decrease; and rapid growth for St > 8. The magnitude of -nrn 2 is_a strong 
decreasing function of 5. This result may be explained as follows. -nin 2 is de- 
termined by two simultaneous effects, a mean deformation effect and a turbulent 
diffusion effect. Let us consider, for instance, a sinusoidal flame surface of constant 
amplitude subjected at t = 0 to a simple shearing motion: — njn 2 first increases, 
saturates at a time that scales with St, and then decreases as the flame elements 
are turned away from the direction of maximum positive strain rate (inclined at 
45 degrees from the mean flow direction) by the rotation component of the applied 
shear. In the simulations, this mean deformation effect is coupled with an increase 
of the turbulent flame thickness due to turbulent diffusion that is constantly coun- 
teracting the flattening of the flame due to the mean rotation. The characteristic 
time scale of the mean deformation effect is 1/5. The characteristic time scale of 
the turbulent diffusion effect is a turbulent time scale that is likely to scale with 
1/v^. The largest magnitude for -nTn 2 is obtained for cases where the turbulence 
dominates the mean deformation, that is for small values of 5. 

Hence, in (25) the increase of dUi/dx 2 with 5 is cancelled by a corresponding 
decrease of the coefficient /?. Figure 8(b) illustrates this result and shows that the 
resulting strain rate is more than one order of magnitude smaller than 5. Also, 
in all cases, a is always 5 to 10 times larger than /? and both terms in (25) have 
comparable contributions to At- 
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3. Conclusion 

Direct simulations of premixed flames in turbulent shear flow are used in this 
study to examine the effect of a mean shear motion on the flame surface dynamics. 
The shear is uniform in the unburnt gas, and simulations are performed for different 
values of the mean shear rate, 5. A detailed analysis of the rate of change of flame 
surface densities is performed, based on the exact E-equation. In particular, the 
time evolution of the total flame stretch, k, is compared for different values of S. 
In the expression for the total flame stretch, the production term_by hydrodynamic 
straming is decomposed into two components: a component, A r , that depends 
explicitly on the mean rate of deformation; and a component, a?, that is determined 
by the turbulent flow field. 

The data show that while k scales with the total shear, St, it does not depend 
on the turbulence-to-mean shear time scale ratio, Sr 0 . This is a surprising result 
since the turbulence is strongly influenced by this parameter. This result implies 
that the principal effect of changing the mean shear rate, 5, is to re-scale the time 
evolution of flame stretch. For instance, the results show that for 0 < St < 10, the 
rate of change of flame surface area is small; in contrast, for St > 10, a transition 
to a phase of rapid growth is observed. The main effect of 5 is then to make that 
transition occur faster. Note also that, while flame stretch scales with St, since 
the total flame surface area, {S v ), is related to the time integral of stretch, (Sv) is 
found to scale both with St and 5. 

The relative insensitivity of k to changes in Sr 0 was further studied by examining 
the direct contribution of the applied shear to A T . This contribution is always small 
in the simulations. It is shown to scale with the product of the magnitude of the 
mean velocity gradient and a coefficient that depends on the flame wrinkling. The 
magnitude of that coefficient is determined by competing effects: a positive turbu- 
ent effect and an initially positive and subsequently negative effect controlled by 
the mean rotation component of the applied shear. It is argued that this coefficient 
scales with the inverse of a characteristic turbulent time divided by S so that at 
lgh shear rates, the increased mean velocity gradient is multiplied by a decreased 
coefficient and there is no net variation of At* 

Finally, the question of which term in the E-equation is responsible for the tran- 
sitionto rapid growth of flame surface area, at St > 10, is studied by comparing 5? 
and At. Results show that a T is increasingly dominant in the expression for flame 
stretch; the main effect of the applied shear is to increase the turbulent straining 
motions, and A T may be neglected in the models. 
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